flutter of darrieus wind turbine blades 


Norman D. Ham 

Department of Aeronautics and Astronautics 
Massachusetts Institute of Technology 
Cambridge/ Mass. 02139 

ABSTRACT 


The testing of Darrieus wind turbines has indicated that under certain 
conditions, serious vibrations of the blades can occur, involving flatwise 
bending, torsion, and chordwise bending. It is the purpose of this paper 
to develop a theorr ‘:ical method of predicting the aeroelastic stability of 
the coupled bending and torsional motion of such blades with a view to 
determining the cause of these vibrations and a means of suppressing them. 

I. Introduction 

The troposkien type of vertical-axis wind turbine. Reference 1, embodies 
flexible curved blades connected at each end to the top and bottom of a verti- 
cal shaft which permits the blades to rotate about the vertical axis under the 
action of the wind. The curve of each blade is ideally that of a troposkien? 
i.e., the shape taken by a flexible cable of uniform density and cross section 
whose ends are attached to two points on a vertical axis, when it is spun at 
a constant angular velocity about the vertical axis. In this manner the bend- 
ing stresses in the blades are eliminated. 

The testing of such devices has indicated that under certain conditions, 
serious vibrations of the blades can occur, involving flatwise bending, torsion, 
and chordwise bending. It is the purpose of this report to develop a theoret- 
ical method of predicting the aeroelastic stability of the coupled bending and 
torsional motion of such blades with a view to determining the cause of these 
vibrations and a means of suppressing them. The present analysis is an exten- 
sion of that of Reference 2. 

Flatwise Bending Vibration of the Ideal Troposkien 

Consider the rotating troposkien- type rotor blade shown in Figure 1. 

The bending moment at x is given by 

M(x) = ^ {[y^+y(5)in^ - y(C)} m(C) (x-C) - H(y(o) - y(x)] 

where 

Q = rotor rotational speed 

m(5) == rotor blade mass at ^ per unit of x-distance 

H = blade tension force at mid-point between blade ends 

Differentiating the bending moment equation twice with respect to x, 

2^ ~ [y +y(x)] - y(x) m(x) + H ^ 

° 3x^ 
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From the theory of elasticity for the bending of slender beams. 


d y 

M(x) = El — ^ 

3x 


Therefore 

3x^ 3x^ 3x^ 


Then the blade bending equation of motion is 


li- [El -^J - H 1-^ - n^m(x) (y +y) + m(x)y = 0 
3x^ 3x^ 

Note that m(x) dx = m(s) ds, where s = arc length along the blade, and 
m(s) = m, a constant. Then 


m(x) 



m 




2 1/2 
1 


Equation (1) is the equation of the vibrating blade: the steady-state 
case is described and solved in Reference 1. For present purposes, consider 
the equation which results when it is assumed that y = y (9 ) : 


li. [El ^ (^)^ (|^)^ -H ^ (|^)^ - f^^m(9) [y +y(9)] + m(e)y(0) = 0 
302 30 ^ 30 ^ ( 2 ) 

Let y = J, y„(e) g„(t) 


Substituting into Equation (2), there results 
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H ? 
n=l 

^ ^n 
302 

30 ^ 

^n 

2 oo 

Q m [y + I, y g ] 
o n=l n ^n 

+ m ?, y„ g„ = 0 
n=l n n 

(3) 
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Assume free viljration of the blade in the nth flatwise mode in a vacuum: 


n n 


ico t 
n 


.. 2* n 2 

g = -0) g e = -cj g 

n n^’n n^n 


Substituting into the nth term of Equation (3), i.e., the equation for 
free vibrations in the nth mode, 

.2 sV .n 2 2 .p 2 ^ 

2 2~* ^ ^ ^ ^ 9 - S^m[y+y g]=mu) y g 

36^ 9 q 2 ax n dx 3 x ^n '^o ^n ^n n ^n ^n 

(4) 

Then including all modes, substitute Equation (4) into Equation (3) to obtain 


S 00 2 

m 9 ^ + m Z 03 y g =0 

n=l n n n=l n -'n n 


(5) 


This is the equation for flatwise free vibration in a vacuum of the rotating 
ideal troposkien blade. 

Flatwise Bending Vibration of the Approximate Troposkien 

Reference 3 describes an approximation to the troposkien shape consisting 

of straight and circular- arc segments. In this instance, the origin of the axis 

system shown in Figure 1 is located at the center of the circular-arc segment, 

and the distance y of the origin to the axis of rotation is determined by the 
o 

geometry of the straight segments. 

Reference 4 suggests the following approximate mode shapes for the flatwise 
vibration of non-rotating circular arcs of radius R with pinned ends: 


_ • /ll'l \ r\ 

y = R sin ( — )9 
■^n a 


(6) 


This expression was shown in Reference 2 to be an approximate solution to 
the eigenvalue equation describing the vibration of rotating circular arcs, for 
the case a = tt and y =0. 
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The vibration frecjuencies are given by the expression 


0 9 2 

0) ss (jj + K 

n Hq n 


( 7 ) 


where a) and K must be determined experimentally, or numerically by finite- 

o _ 2 

element analysis. For the special case of Reference 2, - n - 1. 

Then Equations (5), (6), (7) describe approximately the flatwxse free 

bending vibration in a vacuum of the circular-arc segment of the rotating 
approximate troposkien blade. 


4. Torsion/Chordwise Bending Vibration of the Ap proximate Troposkien 

Reference 5 considers the torsion/chordwise bending vibration of non-rotating 
circular arcs. The forces due to rotation are shown in Figure 2. Adding these 
forces to the equations of Reference 5, in the nomenclature of the present analysis 
the torsion/chordwise bending equations become, for the circular-arc segment of 
the approximate troposkien blade. 


I V ^ M , X ttiR 

Z - r 6” - k (z" + rB’ ' ) = ^ 


(z - - 2f2y - 


mR 


(8) 


2 * • - r 6 + k(RB* ' + z* ' ) = 0 


(9) 


where BO) = blade torsional displacement 

z(0) = blade chordwise bending displacement 
R = radius of circular arc 
k = GK/EI^ 

GK = torsional stiffness of blade section 

EI^ = chordwise bending stiffness of blade section 

m *= rotor blade mass per unit length of blade 

Note that flatwise and torsion/chordwise bending are coupled dynamically 
only by the Coriolis forces 2Qym and 2fizm. 
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Let 


z 


nil 


z (0) f (t) 
n n 


6 


nil 


6„{0) f (t) 
n n 


Substituting into Equation (8), and neglecting the Coriolis coupling term 
at present, there results 


nil 



R 
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Assume free vibration of the blade in the nth torsion/chordwise bending mode in 
a vacuum: 
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Substituting into the nth term of Equation (10), i.e., the equation for free 
vibrations in the nth mode. 



IV 

n 


f 

n 


Rg* * f 
n n 


k(z* * 
n 


f + Rg* * 
n n 


f ) ] + z f + ^ z " 

n n n ^2 n 


f 

n 


= -m z (i)“ f 
n n n 


( 11 ) 


Then including all modes, substitute Equation (11) into Equation (10) to obtain 


m 



z f + m co^ z f 
n n n=l n n n 


0 


(12) 


This IS the equation for torsion/chordwise free vibration in a vacuum of the 
rotating approximate troposkien blade, not including Coriolis coupling with 
flatwise vibration of the blade. 
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Reference 2 suggests the following approximate mode shapes for the 

torsion/chordwise bending vibration of non-rotating circular arcs for the 

case a = TT and y =0; 

o 

z = R sin n0 
n 

6 = sin n0 

n 

Cn 

These expressions were shown in Reference 2 to be approximate solutions to the 
eigenvalue equation describing the vibration of rotating circular arcs. Note 
that they satisfy only the pinned-end boundary condition, and are therefore 
even more approximate for other end conditions, such as the elastic restraint 
of the curved segment by the straight segments of the approximate troposkien 
blade. 

For the case a 7T, the following approximate mode shapes will be used, 
as for the flatwise modes: 

z = R sin (— )0 (13) 

n a 

B = sin (— )0 (i4) 

n a 

C, 

The vibration frequencies are given by the expression 

+ ic (15) 

n n n 
o 

where O) and K must be determined experimentally, or niomerically by finite- 
n n 

^ - 2 
element analysis. For the special case of Reference 2, = n -1. 

Then Equations (12), (13), and (14) describe approximately the torsion/ 

chordwise bending free vibration in a vacuum of the circular-arc segment of the 

rotating approximate troposkien blade. 
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Note that Equation (9) expresses the structural coupling between torsion 
and chordwise bending of the circular arc. Substitution of Equations (13) and 
(14) into Equation (9) yields the nth mode structural coupling coefficient for 
circular arcs having pinned ends. 


C 

n 


z 

n 




1+k 


l+k(iHL)2 

a 


( 16 ) 


This approximate relation was investigated experimentally in Reference 6 and 
found to be valid for the troposkien blade with pinned ends for the case n = 2, 
if an effective value = 2.8 is used to account for the flexibility of the 

straight segments of the blade. 

For the case of clamped ends, the approximate nth mode structural coupling 
coefficient for circular arcs with clamped ends is shown in Reference 6 to be 


C = 
n 


0.75 


. 2 ,n7T Q 

sin (— ) 0 


.25] (^) 


Itk 

a 


where the effective value =2.8 again applies. 

Equation (16) is believed to be a reasonable approximation for the structural 
coupling coefficient of practical blade configurations in the absence of a more 
precise determination by finite element analysis of the blade. 

5, Blade Flutter Equations 

Adding the Coriolis forces due to blade bending, and the blade aerodynamic 
llff force, as shown in Figures 2 and 3, Equations (5) and (12) become 


^ y_ 9 + X mii) V q + 2mi^ z f = ttt sino 

n~l n n n=i n “n 'n n=i n a n dO 


(17) 


■ t 2mn y g + Z, mz f + ^mca^z f = 0 
n~i n n n=l n n n=l n n n 


(18) 
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where chordwise aerodynamic forces are neglected. 

Using Equations (6) and (13), multiplying Equations (17) and (18) by 

R sin(^^)9, and integrating from 0 to ot, the coupled equations of motion in 
Ot 

the nth mode are obtained: 


I g + I o)^g + 21 U f - sin6 sin(^)0 ^ d9 

n n n n ^n n n o a a« 


(19) 


-21 flg+If+IO)f=0 
n n n n n n n 


( 20 ) 


where I = m sin^ (^) 0 d0 
no a 


and noting that ^ sin(^)0 sin(^)0 d9 = 0. 

Following Reference 7 and Figure 3, the aerodynamic lift for the wind 
turbine rotating in still air at constant rotational speed 0, is 

R df “ ■ 8 ■ 2 

where p = air density 

a = blade section lift curve slope 
c *= blade chord 

B *= blade twist angle (postive nose down) (see Fig. 2) 

= outer blade undisturbed shape = Yq ^ sin6 
C(k) = Theodorsen's lift deficiency function 

X « distance from section aerodynamic center to blade elastic 

A 

2 ucis, positive when A.C. is forward. 

Neglecting the B terms since c « R, multiplying by R sin(^)0 and dividing 
2 

I n , there results 
n 
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1 /•TT • X ■ / \ ft dli j ^ 

r / sin6 sin(— )0 = 

j ® ^ 

”*■11 


- m.g - m^C f 
g n 3 n n 


R6 

where C = — - as before. 
n z 


Then 


in. = ^ C(k) 

g 2 


/^(y /R + sin 6 ) sin ^5 sin ^(— )9 d 0 
o Q O 

^sin ^(^)0 d 0 


^ Y 


^ (y /R + sin 6 )^ sin 6 sin ^(— )0 d 0 


mg = C(k) 


f^sin^ (— ) 0 d0 

o a 


where y = 


C(k) = typical value of C (k) evaluated at y^ = y^ + R sin6 (see Appendix), 
and the contributions of modes other than the nth are neglected. 


Dividing Equations (23) and (24) by and including the aerodynamic 


terms. 


9n 2 ^ 

— +m. g + 2 t^ + m-C f = 0 

g U n n Q 6 nn 


( 21 ) 


g f ^ 

- 2 ^ + V f = Q 

*2 W* n n 


( 22 ) 


CD 

where = [ (-^) + K^] 


1/2 


- 0 ) 2 1/2 

^n-sf ■ I'-f' * V 

These are the flutter equations for the wind turbine rotating in still air. 
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5. 


Calculation of Flutter Boundaries 


Assume the following solutions to Equations (21) and (22) ; 


g 


n 



vQt 

e 


( 23 ) 


f 

n 


- vQt 
f e 
n 


(24) 


where ^ + i Assume at present that C(k) =1. 

Substitution of Equations (23) and (24) in Equations (21) and (22) 
results in two coupled homogeneous algebraic equations. For a solution to 
exist, the determinant of coefficients must vanish; i.e.. 


(V^ + m.V + V^) (2V + m^C^) 

g n 3 ri 

. 2 - 2 , 
- 2V [V + V ] 

n 


= 0 


Expanding the determinant yields the characteristic equation of the system; 

4 3 2 

AV + BV + CV + DV + E = 0 


where the coefficients are combinations of the system parameters. Then by 
Routh's criteria for system stability 

A, B, C, D, E > 0 (25) 

2 2 

and BCD - AD - B E > 0 (26) 

It is now possible to determine the values of required to satisfy conditions 
(25) and (26) for a given blade configuration. Since 


V 


2 


n 


0) 

(^) 


+ K 

n 
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and 


v2 

n 


0 ) 2 
n 


+ K 

n 


the values of rotational speed at which flutter or divergence will occur {Q^) 

are then specified for given non-rotating flatwise and torsional frequencies 

0) and 0) 
n n 

o o 

In the present case, the characteristic equation is 

^ ^ ?"22 ^2 2^2 
+ m.V'^ + (4 + + V )V + (m.V + C )V + V V =0 

g nn gnpn nn 

Then conditions (25) and (26) determine the rotational speed at flutter, 

Applying condition (25), flutter occurs when 



n 


2 



m. n 

g 


Applying condition (26) , flutter occurs when 
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4mg 



m. 
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(27) 


(28) 


6. Application to an Existing Blade 

The theory developed above is now applied to the preliminary prediction 
of the flutter of the blades of the Sandia seventeen-meter wind turbine. For 
this turbine the parameters for the geometry of the circular arc portion of the 
blade are 


y^ = 9.10 ft. 

R = 18.33 ft. 
a = 2 rad. 

For this case it can be shown that m^ 1.575 m. .Also from 

p g 

Equation (16) for n = 2, k = .0512, =2.8 (Ref. 6): 

C - 4.2 
n 
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Applying Equation (27) , flutter occurs when 

= 13.2 
n 

From Reference (8), for the seventeen-meter blade, 





7.70 


Equating and solving the two expressions for 

^ = 138 rpm 


Applying Equation (28) flutter occurs when 

= 13.2 + 0.768 


From Reference (8) , for the seventeen-meter blade. 


Then 




= 18.4 + 0.768 

2 iip 


Equating and solving the two expressions for v^, 

= 80.4 rpm 


Here Equation (28) is evidently the critical case. 

Final prediction of the flutter speed requires accurate determination of 
the structural coupling coefficient C by finite element analysis. 


7. Conclusions 

(1) Coriolis forces provide the dominant coupling between flatwise bending 
and torsion/chordwise bending of the blade. (Also see Reference 2) . 

(2) Blade flatwise bending velocity and blade torsional displacement 
contribute the dominant aerodynamic forces acting on the blade. 


88 



(3) Accurate determination of the structural coupling between blade 
torsion and blade chordwise bending is essential to blade flutter 
prediction. 

(4) The following parameters determine the rotational speed at flutter 
for troposkien-type blades: 

(a) blade geometry y , R, a (see Figure 1). 

o 

(b) blade torsion/chordwise bending structural coupling coefficient 

for the second mode. 

(c) blade rotating and non-rotating flatwise bending frequencies 

(jL)^ and for the second mode, 
o 

(d) blade rotating and non-rotating torsion/chordwise bending 

frequencies Oi^ and 0)^ for the second mode, 
o 
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DISCUSSION 

Comment: I would like to reemphasize the significance of Professor Ham's 

conclusions that a) quarter-chord mass balancing of the turbine blades 
is unimportant with regard to flutter speed location, and b) increasing 
torsional stiffness is an effective method of increasing flutter speed. 

The impact of the first conclusion on design and fabrication of low cost 
blades is very favorable, and knowledge of the second conclusion is 
valuable to the designer. At the same time, I would like to ^-emphasize 
the threat of encountering flutter in any particular blade design. Based 
on Professor Ham’s analysis, blade property evaluation, field experience, 
and consultation with fabricators, it is relatively easy to design a 
blade that will not flutter . 

Q: Do you agree with R. Reuter that it is relatively easy to design 
Darrieus blades which do not flutter in their operating regime? 

A: Yes. 
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Q: Am I correct in presuming that the forced response to periodic stalling 
should also be independent of the chordwise c.g. of the blade sections, 
and also of the chordwise shear center? 

A: Yes. 

Q: What were your assumptions relative to the aerodynamics; was the 
free stream velocity zero? If not what of the resultant periodic 
coefficient? 

A: In the analysis the free stream velocity was assumed to be zero. If 
this were not so, harmonic airloads would be generated and would 
cause a forced blade response at various integral multiples of rotor 
speed. However, in a linear analysis, this forced response is uncoupled 
with any self-excited motion such as flutter. Experimentally, model 
blade flutter occurred at the same rotor rotational speed whether or 
not there was an incident wind. 
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